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CN . Abstract 

^ I We study properties of programs with monotone and convex constraints. We extend to these 

i formaHsms concepts and results from normal logic programming. They include the notions of 

strong and uniform equivalence with their characterizations, tight programs and Pages Lemma, 
i program completion and loop formulas. Our results provide an abstract account of properties of 

some recent extensions of logic programming with aggregates, especially the formalism of Iparse 
programs. They imply a method to compute stable models of Iparse programs by means of off-the- 
I shelf solvers of pseudo-boolean constraints, which is often much faster than the smodels system. 
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1. Introduction 



\ We study programs with monotone constraints (Marek & Truszczynski, 2004; Marek, Niemela, 

^ ' & Truszczynski, 2004, 2006) and introduce a related class of programs with convex constraints. 

Q . These formalisms allow constraints to appear in the heads of program rules, which sets them apart 

O \ from other recent proposals for integrating constraints into logic programs (Pelov, 2004; Pelov, 

Denecker, & Bruynooghe, 2004, 2006; DeU'Armi, Faber, lelpa, Leone, & Pfeifer, 2003; Faber, 
Leone, & Pfeifer, 2004), and makes them suitable as an abstract basis for formalisms such as Iparse 
programs (Simons, Niemela, & Soininen, 2002). 

We show that several results from normal logic programming generalize to programs with mono- 
^ . tone constraints. We also discuss how these techniques and results can be extended further to the 

H ' setting of programs with convex constraints. We then apply some of our general results to design 

and implement a method to compute stable models of Iparse programs and show that it is often 
much more effective than smodels (Simons et al., 2002). 

Normal logic programming with the semantics of stable models is an effective knowledge repre- 
sentation formalism, mostly due to its ability to express default assumptions (Baral, 2003; Gelfond 
& Leone, 2002). However, modeling numeric constraints on sets in normal logic programming is 
cumbersome, requires auxiliary atoms and leads to large programs hard to process efficiently. Since 
such constraints, often called aggregates, are ubiquitous, researchers proposed extensions of normal 
logic programming with explicit means to express aggregates, and generalized the stable-model se- 
mantics to the extended settings. 

Aggregates imposing bounds on weights of sets of atoms and literals, called weight constraints, 
are especially common in practical applications and are included in all recent extensions of logic 
programs with aggregates. Typically, these extensions do not allow aggregates to appear in the 
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heads of rules. A notable exception is the formalism of programs with weight constraints (Niemela, 
Simons, & Soininen, 1999; Simons et al., 2002), which we refer to as Iparse programs^. 

Lparse programs are logic programs whose rules have weight constraints in their heads and 
whose bodies are conjunctions of weight constraints. Normal logic programs can be viewed as a 
subclass of Iparse programs and the semantics of Iparse programs generalizes the stable-model 
semantics of normal logic programs (Gelfond & Lifschitz, 1988). Lparse programs are one of the 
most commonly used extensions of logic programming with weight constraints. 

Since rules in Iparse programs may have weight constraints as their heads, the concept of one- 
step provability is nondeterministic, which hides direct parallels between Iparse and normal logic 
programs. An explicit connection emerged when Marek and Truszczynski (2004) and Marek et al. 
(2004, 2006) introduced logic programs with monotone constraints. These programs allow aggre- 
gates in the heads of rules and support nondeterministic computations. Marek and Truszczynski 
(2004) and Marek et al. (2004, 2006) proposed a generalization of the van Emden-Kowalski one- 
step provability operator to account for that nondeterminism, defined supported and stable models 
for programs with monotone constraints that mirror their normal logic programming counterparts, 
and showed encodings of smodels programs as programs with monotone constraints. 

In this paper, we continue investigations of programs with monotone constraints. We show that 
the notions of uniform and strong equivalence of programs (Lifschitz, Pearce, & Valverde, 2001; 
Lin, 2002; Turner, 2003; Eiter & Fink, 2003) extend to programs with monotone constraints, and 
that their characterizations (Turner, 2003; Eiter & Fink, 2003) generalize, too. 

We adapt to programs with monotone constraints the notion of a tight program (Erdem & Lifs- 
chitz, 2003) and generalize Fages Lemma (Fages, 1994). 

We introduce extensions of propositional logic with monotone constraints. We define the com- 
pletion of a monotone-constraint program with respect to this logic, and generaUze the notion of a 
loop formula. We then prove the loop-formula characterization of stable models of programs with 
monotone constraints, extending to the setting of monotone-constraint programs results obtained 
for normal logic programs by Clark (1978) and Lin and Zhao (2002). 

Programs with monotone constraints make explicit references to the default negation operator. 
We show that by allowing a more general class of constraints, called convex, default negation can be 
eliminated from the language. We argue that all results in our paper extend to programs with convex 
constraints. 

Our paper shows that programs with monotone and convex constraints have a rich theory that 
closely follows that of normal logic programming. It impUes that programs with monotone and con- 
vex constraints form an abstract generalization of extensions of normal logic programs. In particu- 
lar, all results we obtain in the abstract setting of programs with monotone and convex constraints 
specialize to Iparse programs and, in most cases, yield results that are new. 

These results have practical implications. The properties of the program completion and loop 
formulas, when specialized to the class of Iparse programs, yield a method to compute stable models 
of Iparse programs by means of solvers of pseudo-boolean constraints, developed by the proposi- 
tional satisfiability and integer programming communities (Een & Sorensson, 2003; Aloul, Ramani, 
Markov, & Sakallah, 2002; Walser, 1997; Manquinho & Roussel, 2005; Liu & Truszczynski, 2003). 
We describe this method in detail and present experimental results on its performance. The results 
show that our method on problems we used for testing typically outperforms smodels. 

1. Aggregates in the heads of rules have also been studied recently by Son and Pontelli (2006) and Son, PontelU, and 
Tu (2006). 
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2. Preliminaries 

We consider the propositional case only and assume a fixed set At of prepositional atoms. It does 
not lead to loss of generality, as it is common to interpret programs with variables in terms of their 
propositional groundings. 

The definitions and results we present in this section come from papers by Marek and Trusz- 
czynski (2004) and Marek et al. (2006). Some of them are more general as we allow constraints 
with infinite domains and programs with inconsistent constraints in the heads. 
Constraints. A constraint is an expression A = {X, C), where X C At and C C ViX) {V{X) 
denotes the powerset of X). We call the set X the domain of the constraint A = (X, C) and denote it 
by Dom{A). Informally speaking, a constraint (X, C) describes a property of subsets of its domain, 
with C consisting precisely of these subsets of X that satisfy the constraint (have property) C. 

In the paper, we identify truth assignments (interpretations) with the sets of atoms they assign 
the truth value true. That is, given an interpretation M C ^i, we have M ^ a if and only if a G M. 
We say that an interpretation M C At satisfies a constraint A = {X, C) {M \= A),\i M r^ X e C. 
Otherwise, M does not satisfy A, (M y= A). 

A constraint A = {X, C) is consistent if there is M such that M \= A. Clearly, a constraint 
A = (X, C) is consistent if and only if C 7^ 0. 

We note that propositional atoms can be regarded as constraints. Let a ^ At and M C At. We 
define C{a) = ({a}, {{a}})- It is evident that M \= C{a) if and only if M |= a. Therefore, in the 
paper we often write a as a shorthand for the constraint C{a). 

Constraint programs. Constraints are building blocks of rules and programs. Marek and Trusz- 
czynski (2004) defined constraint programs as sets of constraint rules 

A*^ Al,..., ^fe, not(^fe+i), . . . , not(yl^) (1) 

where A, Ai, . . . , An are constraints and not is the default negation operator. 

In the context of constraint programs, we refer to constraints and negated constraints as literals. 
Given a rule r of the form (1), the constraint (literal) A is the head of r and the set {^1, . . . , 
Aj^, ... , not(Ajfc+i), . . . , not{Ajn)} of literals is the body of r^. We denote the head and the body 
of r by hd{r) and bd{r), respectively. We define the the headset of r, written hset(r), as the domain 
of the head of r. That is, hset{r) = Dom{hd{r)). 

For a constraint program P, we denote by At{P) the set of atoms that appear in the domains of 
constraints in P. We define the headset of P, written hset{P), as the union of the headsets of all 
rules in P. 

Models. The concept of satisfiability extends in a standard way to literals not(^) (M |= not(^) 
if M ^ A), to sets (conjunctions) of literals and, finally, to constraint programs. 
M-applicable rules. Let M C Athe an interpretation. A rule (1) is M-applicable if M satisfies 
every literal in bd{r). We denote by P{M) the set of all M-applicable rules in P. 
Supported models. Supportedness is a property of models. Intuitively, every atom a in a supported 
model must have "reasons" for being "in". Such reasons are M-applicable rules whose heads con- 
tain a in their domains. Formally, let P be a constraint program and M a subset of At{P). A model 
M of P supported if M C hset{P{M)). 

Examples. We illustrate the concept with examples. Let P be the constraint program that consists 
of the following two rules: 

2. Sometimes we view the body of a rule as tlie conjunction of its literals. 
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({c, d, e}, {{c}, {d}, {e}, {c, d, e}}) ^ 

({a, 6}, {{a}, {6}}) ^ {{c,d},{{c},{c,d}}), not {i{e},{{e}})) 

A set M = {a, c} is a model of P as M satisfies the heads of the two rules. Both rules in P are 
M-applicable. The first of them provides the support for c, the second one — for a. Thus, M is a 
supported model. 

A set M' = {a, c, d, e} is also a model of P. However, a has no support in P. Indeed, a only 
appears in the headset of the second rule. This rule is not M'-apphcable and so, it does not support 
a. Therefore, M' is not a supported model of P. A 
Nondeterministic one-step provability. Let P be a constraint program and M a set of atoms. A set 
M' is nondeterministically one-step provable from M by means of P, if M' C hset{P{M)) and 
M' ^ hd(r), for every rule r in P(M). 

The nondeterministic one-step provability operator Tp^ for a program P is an operator on 
V{At) such that for every M C At, Tp'^{M) consists of all sets that are nondeterministically 
one-step provable from M by means of P. 

The operator Pp"' is nondeterministic as it assigns to each M C At a family of subsets of At, 
each being a possible outcome of applying P to M. In general, Tp'^ is partial, since there may be 
sets M such that Tp'^{M) = (no set can be derived from M by means of P). For instance, if 
P{M) contains a rule r such that hd{r) is inconsistent, then Tp'^{M) = 0. 

Monotone constraints. A constraint (X, C) is monotone if C is closed under superset, that is, for 
every VF,y C X, if G C and C y then F G C. 

Cardinality and weight constraints provide examples of monotone constraints. Let X be n finite 
set and let Ck{X) = {Y : Y (1 X, k < \Y\}, where /c is a non-negative integer. Then {X, Ck{X)) 
is a constraint expressing the property that a subset of X has at least k elements. We call it a lower- 
bound cardinality constraint on X and denote it by kX. 

A more general class of constraints are weight constraints. Let X be a finite set, say X = 
{xi, . . . , Xn}, and let w, wi, . . . , be non-negative reals. We interpret each Wi as the weight as- 
signed to Xj. A lower-bound weight constraint is a constraint of the form {X, C^), where con- 
sists of those subsets of X whose total weight (the sum of weights of elements in the subset) is at 
least w. We write it as 

w[xi =Wi,...,Xn= Wn]- 

If all weights are equal to 1 and w is an integer, weight constraints become cardinality con- 
straints. We also note that the constraint C{a) is a cardinality constraint l{a} and also a weight 
constraint l[a = 1]. Finally, we observe that lower-bound cardinality and weight constraints are 
monotone. 

Cardinality and weight constraints (in a somewhat more general form) appear in the language of 
Iparse programs (Simons et al., 2002), which we discuss later in the paper. The notation we adopted 
for these constraints in this paper follows the one proposed by Simons et al. (2002). 

We use cardinality and weight constraints in some of our examples. They are also the focus of 
the last part of the paper, where we use our abstract results to design a new algorithm to compute 
models of Iparse programs. 

Monotone-constraint programs. We call constraint programs built of monotone constraints — 

monotone-constraint programs or programs with monotone constraints. That is, monotone-constraint 
programs consist of rules of rules of the form (1), where A, Ai, . . . , Am are monotone constraints. 
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From now on, unless explicitly stated otherwise, programs we consider are monotone-constraint 
programs. 

2.1 Horn Programs and Bottom-up Computations 

Since we allow constraints with infinite domains and inconsistent constraints in heads of rules, the 
results given in this subsection are more general than their counterparts by Marek and Truszczyiiski 
(2004) and Marek et al. (2004, 2006). Thus, for the sake of completeness, we present them with 
proofs. 

A rule (1) is Horn ifk = m (no occurrences of the negation operator in the body or, equivalently, 
only monotone constraints). A constraint program is Horn if every rule in the program is Horn. 

With a Horn constraint program we associate bottom-up computations, generaUzing the corre- 
sponding notion of a bottom-up computation for a normal Horn program. 

Definition 1. Let P be a Horn program. A P-computation is a (transfinite) sequence {Xa) such 
that 



2. for every ordinal number a, C X^+i and X^+i G Tp (Xq), 

3. for every limit ordinal a, X^ = U/3<a ^p- 

Let t = {Xa) be a P-computation. Since for every (5 < (3', X^ C X^i C At, there is a least 
ordinal number at such that X^t+i = > other words, a least ordinal when the P-computation 
stabilizes. We refer to at as the length of the P-computation t. 

Examples. Here is a simple example showing that some programs have computations of length 
exceeding u) and so, the transfinite induction in the definition cannot be avoided. Let P be the 
program consisting of the following rules: 



where Xi = {ao, . . . aj}, < i, and X^ = {ao, ai, . . .}. Since the body of the last rule con- 
tains a constraint with an infinite domain X^o, it does not become appUcable in any finite step of 
computation. However, it does become appUcable in the step u and so, a G -'^w-t-i- Consequently, 



For a P-computation t = {X^), we call (J^ X^ the result of the computation and denote it by 
Rt. Directly from the definitions, it follows that Rt = Xat- 

Proposition 1. Let P be a Horn constraint program and t a P-computation. Then Rt is a supported 
model of P. 

Proof. Let M = Rthe the result of a P-computation t = {Xa). We need to show that: (1) M is a 
model of P; and (2) M C hset{P{M)). 

(1) Let us consider a rule r & P such that M \= bd{r). Since M = Rt = X^t (where at is the 
length of t), Xat N bd{r). Thus, Xat+i \= hd{r). Since M = X^^+i, M is a model of r and, 
consequently, of P, as well. 



1. Xo = 0, 




A 
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(2) We will prove by induction that, for every set in the computation t, X^ C hset{P{M)). The 
base case holds since = C hset{P{M)). 

If a = /? + 1, then X^ G T^'^iX/s). It follows that X^ C /isei(P(X^)). Since P is a Horn 
program and X/s C M, hset{P{Xp)) C hset{P{M)). Therefore, C hset{P{M)). 

If a is a limit ordinal, then Xa = U^<a -^iS- induction hypothesis, for every /3 < a, 

Xf) C hset{P{M)). Thus, C hset{P{M)). By induction, M C hset{P{M)). □ 

Derivable models. We use computations to define derivable models of Horn constraint programs. 

A set M of atoms is a derivable model of a Horn constraint program P if for some P-computation 
t, we have M = Rt- By Proposition 1, derivable models of P are supported models of P and so, 
also models of P. 

Derivable models are similar to the least model of a normal Horn program in that both can be 
derived from a program by means of a bottom-up computation. However, due to the nondeterminism 
of bottom-up computations of Horn constraint programs, derivable models are not in general unique 
nor minimal. 

Examples. For example, let P be the following Horn constraint program: 

P={l{a,b} ^} 

Then {a}, {b} and {a, b} are its derivable models. The derivable models {a} and {b} are minimal 
models of P. The third derivable model, {a, b}, is not a minimal model of P. A 

Since inconsistent monotone constraints may appear in the heads of Horn rules, there are Horn 
programs P and sets X C At, such that Tp^{X) = 0. Thus, some Horn constraint programs have 
no computations and no derivable models. However, if a Horn constraint program has models, the 
existence of computations and derivable models is guaranteed. 

To see this, let M be a model of a Horn constraint program P. We define a canonical compu- 
tation t^'^ = {Xa'^) by specifying the choice of the next set in the computation in part (2) of 
Definition 1. Namely, for every ordinal P, we set 

X^:^ = hset{P{X^'''))nM. 

That is, we include in Xa'^ all those atoms occurring in the heads of xj'^-applicable rules that 
belong to M. We denote the result of t^'^ by Can{P,M). Canonical computations are indeed 
P-computations. 

Proposition 2. Let P be a Horn constraint program. If M C. At is a model of P, the sequence t^'^ 
is a P -computation. 

Proof. As P and M are fixed, to simplify the notation in the proof we will write Xa instead of 

To prove the assertion, it suffices to show that (1) hset{P{Xa)) n M e T]^'^{Xa), and (2) 
Xa C hset{P{Xa)) n M, for every ordinal a. 

(I) Let X C M and r G P{X). Since all constraints in bd{r) are monotone, and X \= bd{r), M |= 
bd{r), as well. From the fact that M is a model of P it follows now that M \= hd{r). Consequently, 
M n hset{P{X)) \= hd{r) for every r G P{X). Since M n hset{P{X)) C hset{P{X)), 

M n hset{P{X)) G TP'^(X). 
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Directly from the definition of the canonical computation for P and M we obtain that for every 
ordinal a, ^ M. Thus, (1), follows. 

(2) We proceed by induction. The basis is evident as Xq = 0. Let us consider an ordinal a > 
and let us assume that (2) holds for every ordinal P < a.\i a = (3 + 1, then Xa = Xj^j^i = 
hset{P{Xp)) n M. Thus, by the induction hypothesis, Xp C X^. Since P is a Horn constraint 
program, it follows that P{Xp) C P{Xa). Thus 

X^ = X/s+i = hset{P{Xfi)) n M C hset{P{Xa)) n M. 

If a is a hmit ordinal then for every (3 < a, X^ C X^ and, as before, also P{Xp) C P[Xa). Thus, 
by the induction hypothesis for every (5 < a, 

Xp C hset{P{X/))) n M C hset{P{Xa)) n M, 

which implies that 

Xa=\JXpC hset{P{Xa)) n M. 

□ 

Canonical computations have the following ^/join? property. 
Proposition 3. Let P be a Horn constraint program. For every model M ofP, we have 

hset{P{Can{P, M))) n M = Can{P, M). 

Proof. Let a be the length of the canonical computation t^'^ . Then, X^^ = X^'^ = Can{P, M). 
Since X^+i = hset{Xa) H M, the assertion follows. □ 

We now gather properties of derivable models that extend properties of the least model of normal 
Horn logic programs. 

Proposition 4. Let P be a Horn constraint program. Then: 

L For every model M of P, Can(P, M) is a greatest derivable model of P contained in M 

2. A model M of P is a derivable model if and only ifM = Can{P, M) 

3. If M is a minimal model ofP then M is a derivable model of P. 

Proof (1) Let M' be a derivable model of P such that M' C M. Let T = (Xa) be a P-derivation 
such that M' = Rt. We will prove that for every ordinal a, X^ C Xa'^ ■ We proceed by transfinite 
induction. Since Xq = Xq''^ = 0, the basis for the induction is evident. Let us consider an ordinal 
a > and assume that for every ordinal (3 < a, Xf^ C X^'^ . 

If a = /5 + 1, then Xa G T^'^{Xp) and so, Xa C hset{P{Xp)). By the induction hypothesis 
and by the monotonicity of the constraints in the bodies of rules in P, Xa ^ hset{P{X^^ '^^ ) )- Thus, 
since C P^ = M' C M, 

Xa C hset{P{Xp'')) n M = X^fl = X^'^. 
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The case when a is a Umit ordinal is straightforward as Xa = U/3<a ^^^^ ^a'^ = U/3<a -^^'^ ■ 

(2) (<s=) If M = Can{P, M), then M is the result of the canonical P-derivation for P and M. In 
particular, M is a derivable model of P. 

(^) if M is a derivable model of P, then M is also a model of P. From (1) it follows that 
Can{P, M) is the greatest derivable model of P contained in M. Since M itself is derivable, 

M = Can{P,M). 

(3) From (1) it follows that Can{P, M) is a derivable model of P and that Can{P, M) C M. Since 
M is a minimal model, Can{P, M) = M and, by (2), M is a derivable model of P. □ 

2.2 Stable Models 

In this section, we will recall and adapt to our setting the definition of stable models proposed and 
studied by Marek and Truszczynski (2004) and Marek et al. (2004, 2006) Let P be a monotone- 
constraint program and M a subset of At(P). The redact of P, denoted by P*^, is a program 
obtained from P by: 

1. removing from P all rules whose body contains a literal Yvot{B) such that M \= B; 

2. removing Uterals not(P) for the bodies of the remaining rules. 

The reduct of a monotone-constraint program is Horn since it contains no occurrences of default 
negation. Therefore, the following definition is sound. 

Definition 2. Let P be a monotone-constraint program. A set of atoms M is a stable model ofP if 
M is a derivable model of P^ . We denote the set of stable models of P by St{P). 

The definitions of the reduct and stable models follow and generalize those proposed for normal 
logic programs, since in the setting of Horn constraint programs, derivable models play the role of 
a least model. 

As in normal logic programming and its standard extensions, stable models of monotone- 
constraint programs are supported models and, consequently, models. 

Proposition 5. Let P be a monotone-constraint program. If M C. At{P) is a stable model of P, 
then M is a supported model of P. 

Proof. Let M be a stable model of P. Then, M is a derivable model of P^ and, by Proposition 1, 
M is a supported model of P^ . It follows that M is a model of P^ . Directly from the definition 
of the reduct it follows that M is a model of P. 

It also follows that M C hset{P^{M)). For every rule r in P^^(M), there is a rule r' in P(M), 
which has the same head and the same non-negated literals in the body as r. Thus, /tsei (P^(M)) C 
hset{P{M)) and, consequently, M C hset{P{M)). It follows that M is a supported model of 
P. □ 

Examples. Here is an example of stable models of a monotone-constraint program. Let P be a 
monotone-constraint program that contains the following rules: 

2{a,b,c} ^ l{a,d},not{l{c}) 
l{b,c,d} ^ l{a},not{3{a,b,d})) 
l{a}^ 
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Let M = {a, b). Therefore, M ^ l{c} and M ^ 3{a, 6, d}. Hence the reduct contains the 
following three Horn rules: 

2{a, h, c} ^ l{a, d} 
l{&,c,4 ^ l{a} 
l{a}^ 

Since M = {a, 6} is a derivable model of , M is a stable model of P. 

Let M' = {a, b, c}. Then M' \= l{c} and M ^ 3{a, b, d}. Therefore, the reduct P^' contains 
two Horn rules: 

l{b,c,d} ^ l{a} 
l{a}^ 

Since M' = {a, b, c} is a derivable models of P^' , M' is also a stable model of P. We note that 
stable models of a monotone-constraint program, in general, do not form an anti-chain. A 
If a normal logic program is Horn then its least model is its (only) stable model. Here we have 
an analogous situation. 

Proposition 6. Let P be a Horn monotone-constraint program. Then M C At{P) is a derivable 
model ofP if and only ifM is a stable model of P. 

Proof. For every set M of atoms P = P^. Thus, M is a derivable model of P if and only if it is a 
derivable model of P^ or, equivalently, a stable model of P. □ 

In the next four sections of the paper we show that several fundamental results concerning 
normal logic programs extend to the class of monotone-constraint programs. 

3. Strong and Uniform Equivalence of Monotone-constraint Programs 

Strong equivalence and uniform equivalence concern the problem of replacing some rules in a logic 
program with others without changing the overall semantics of the program. More specifically, 
the strong equivalence concerns replacement of rules within arbitrary programs, and the uniform 
equivalence concerns replacements of all non-fact rules. In each case, the stipulation is that the 
resulting program must have the same stable models as the original one. Strong (and uniform) 
equivalence is an important concept due to its potential uses in program rewriting and optimization. 

Strong and uniform equivalence have been studied in the hterature mostly for normal logic 
programs (Lifschitz et al., 2001; Lin, 2002; Turner, 2003; Eiter & Fink, 2003). 

Turner (2003) presented an elegant characterization of strong equivalence of smodels programs, 
and Eiter and Fink (2003) described a similar characterization of uniform equivalence of normal 
and disjunctive logic programs. We show that both characterizations can be adapted to the case of 
monotone-constraint programs. In fact, one can show that under the representations of normal logic 
programs as monotone-constraint programs (Marek et al., 2004, 2006) our definitions and charac- 
terizations of strong and uniform equivalence reduce to those introduced and developed originally 
for normal logic programs. 
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3.1 M-maximal Models 

A key role in our approach is played by models of Horn constraint programs satisfying a certain 
maximality condition. 

Definition 3. Let P be a Horn constraint program and let M be a model of P. A set N C. M such 
that N is a model ofP and Mnhset{P{N)) C N is an M-maximal model ofP, written N \=m P- 

Intuitively, N is an M-maximal model of P if satisfies eacfi rule r € P{N) "maximally" 
with respect to M. That is, for every r G P{N), N contains all atoms in M that belong to hset{r) 
— the domain of the head of r. 

To illustrate this notion, let us consider a Horn constraint program P consisting of a single rule: 

l{p,q,r} ^ l{s,t}. 

Let M = {p, g, s, t} and N = {p, q, s}. One can verify that both M and A'" are models of P. 
Moreover, since the only rule in P is iV-applicable, and M n {p, q, r} C N, N is an M-maximal 
model of P. On the other hand, N' = {p, s} is not M-maximal even though A^' is a model of P and 
it is contained in M. 

There are several similarities between properties of models of normal Horn programs and M- 
maximal models of Horn constraint programs. We state and prove here one of them that turns out to 
be especially relevant to our study of strong and uniform equivalence. 

Proposition 7. Let P he a Horn constraint program and let M be a model of P. Then M is an 

M-maximal model of P and Can{P, M) is the least M -maximal model of P. 

Proof. The first claim follows directly from the definition. To prove the second one, we simplify the 
notation: we will write for Can{P, M) and for Xa'^ . 

We first show that AT is an M-maximal model of P. Clearly, N C M. Moreover, by Proposition 
3, hset{P{N)) nM = N. Thus, A^ is indeed an M-maximal model of P. 

We now show A^ is the least M-maximal model of P. 

Let A^' be any M-maximal model of P. We will show by transfinite induction that A'' C N'. 
Since Xq = 0, the basis for the induction holds. Let us consider an ordinal a > and let us assume 
that X/3 C N', for every /? < a. To show A^ C A^', it is sufficient to show that Xa C A^'. 

Let us assume that a = ^5 + 1 for some /3 < a. Then, since X/^ C N' and P is a Horn constraint 
program, we have P(X^) C P(N'). Consequently, 

Xa = X^+i = hset{P{Xp)) n M C hset{P{N')) n M C AT', 

the last inclusion follows from the fact ihatN' is an M-maximal model of P. 

If a is a limit ordinal, then Xa = U p^a -^P inclusion Xa C N' follows directly from 

the induction hypothesis. □ 

3.2 Strong Equivalence and SE-models 

Monotone-constraint programs P and Q are strongly equivalent, denoted by P =s Q, if for every 
monotone-constraint program R, P U R and Q U R have the same set of stable models. 

To study the strong equivalence of monotone-constraint programs, we generaUze the concept of 
an SE-model due to Turner (2003). 
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There are close connections between strong equivalence of normal logic programs and the logic 
here-and-there. The semantics of the logic here-and-there is given in terms of Kripke models with 
two words which, when rephrased in terms of pairs of interpretations (pairs of sets of propositional 
atoms), give rise to SE-models. 

Definition 4. Let P be a monotone-constraint program and let X, Y be sets of atoms. We say that 
{X,Y) is an SE-model of P if the following conditions hold: (1) X C Y; (2) Y \= P; and (3) 
X \=Y P^. We denote by SE{P) the set of all SE-models of P. 

Examples. To illustrate the notion of an SE-model of a monotone-constraint program, let P consist 
of the following two rules: 

2{p, q, r} ^ l{q, r}, not(3{p, q, r})} 
l{p,s} ^ l{p,r},noti2{p,r}) 

We observe that M = {p, q} is a model of P. Let N = $. Then N C M md P^{N) is empty. It 
follows that M n hset{P^{N)) = C N aad so, N \=m ■ Hence, {N, M) is an SE-models 
of P. 

Next, let N' = {p}. It is clear that A^' C M. Moreover, P^{N') = {l{p,s} ^ l{p,r}}. 
Hence M n hset{P^{N')) = {p} C N' and so, N' \=m P^ ■ That is, {N' , M) is another SE- 
model of P. A 

SE-models yield a simple characterization of strong equivalence of monotone-constraint pro- 
grams. To state and prove it, we need several auxiliary results. 

Lemma 1. Let P be a monotone-constraint program and let M bea model of P. Then (M, M) and 

{Can{P^, M), M) are both SE-models of P. 

Proof. The requirements (1) and (2) of an SE-model hold for (M, M). Furthermore, since M is a 
model of P, M 1= P^ . Finally, we also have hset{P{M)) nM CM. Thus, M \=m P^ ■ 

Similarly, the definition of a canonical computation and Proposition 1, imply the first two re- 
quirements of the definition of SE-models for {Can{P^ , M),M). The third requirement follows 
from Proposition 7. □ 

Lemma 2. Let P and Q be two monotone-constraint programs such that SE{P) = SE{Q). Then 
St{P) = St{Q). 

Proof If M G St{P), then M is a model of P and, by Lemma 1, (M,M) G SE{P). Hence, 
(M, M) e SE{Q) and, in particular, M |= Q. By Lemma 1 again, 

{Can{Q^,M),M) e SE{Q). 

By the assumption, 

{Can{Q^,M),M) e SE{P) 

and so, Can{Q^'^ , M) or, in other terms, Can{Q^, M) is an M-maximal model of P^ . 

Since M e St{P), M = Can{P^,M). By Proposition 7, M is the least M-maximal model 
of P^. Thus, M C Can{Q^,M). On the other hand, we have Can{Q'^,M) C M and so, 
M = Can{Q^, M). It follows that M is a stable model of Q. The other inclusion can be proved 
in the same way. □ 
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Lemma 3. Let P and R be two monotone-constraint programs. Then SE{P U R) = SE{P) fl 
SE{R). 

Proof. The assertion follows from the following two simple observations. First, for every set Y of 
atoms, Y \= {P U R) if and only if y ^ P and Y \= R. Second, for every two sets X and Y of 
atoms, X\=Y (PUR)^ if and only if X \=y P^ and X \=y R^ . □ 

Lemma 4. Let P, Q be two monotone-constraint programs. If P =s Q, then P and Q have the 
same models. 

Proof. Let M be a model of P. By r we denote a constraint rule (M, {M}) ^ . Then, M G 
St{P U {r}). Since P and Q are strongly equivalent, M e St{Q U {r}). It follows that M is a 
model of Q U {r} and so, also a model of Q. The converse inclusion can be proved in the same 
way. □ 

Theorem 1. Let P and Q be monotone-constraint programs. Then P =s Q if and only ifSE{P) = 
SE{Q). 

Proof. (<;=) Let R be an arbitrary monotone-constraint program. Lemma 3 implies that SE{P U 
R) = SE{P) n SE{R) and SE{Q U R) = SE{Q) n SE{R). Since SE{P) = SE{Q), we have 
that SE{P \J R) = SE{Q U R). By Lemma 2, P U R aad Q U R have the same stable models. 
Hence, P =s Q holds. 

{=>) Let us assume SE{P) \ SE{Q) / and let us consider {X, Y) G SE{P) \ SE{Q). It follows 
that X C y and y ^ P. By Lemma 4,Y \=Q. Since {X, Y) ^ SE{Q), X Q^. It follows 
that X ^ or hset{Q^ (X)) n y ^ X. In the first case, there is a rule r G Q^{X) such that 
X ^ hd{r). Since X QY and is a Horn constraint program, r G {Y). Let us recall that 
y [= Q and so, we also have Y ^ Q^. It follows that Y |= /id(r). Since hset(r) C /isei(Q^(X)), 
y n hset{Q^{X)) \= hd{r). Thus, hset{Q^ {X)) r\Y % X (otherwise, by the monotonicity of 
hd{r), we would have X \= hd{r)). 

The same property holds in the second case. Thus, it follows that 

{hset{Q^ {X)) r\Y)\X ^ ^. 

We define 

X' = {hset{Q^{X)) n y) \ X 
Let i? be a constraint program consisting of the following two rules: 

(X, {X}) ^ 
{Y,{Y})^{X\{X'}). 

Let us consider a program Qo = Q U P. Since y ^ Q and X C y, y ^ Qq- Thus, Y ^ Ql and, 
in particular, Can{Q^ ,Y) is well defined. Since R C , X C Can{Q^ ,Y). Thus, we have 

hset{Ql{X)) n y C hset{Ql{Can{Ql, Y))) r\Y = Can{Ql, Y) 

(the last equaUty follows from Proposition 3). We also have Q QQq and so, 

X' c hset{Q^{x)) n y c hset{Ql{x)) n y 
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Thus, X' C Can{QQ , Y). Consequently, by Proposition 3,Y C Can{Q^ , Y). Since Can{Q^, Y) 
CY,Y = Can{Ql, Y) and so, Y € St{Qo). 

Since P and Q are strongly equivalent, Y G St{Po), where Po = P U R. Let us recall that 
{X,Y) G SE{P). By Proposition 7, Can{P^,Y) is a least F-maximal model of P^. Since X 
is a y-maximal model of P (as X \=y P^), it follows that Can{P^ , Y) C X. Since X' g X, 
Can{P^, Y) C X. Finally, since X' CY,Y ^ X. Thus, y 7^ Can{P^,Y), a contradiction. 

It follows that SE{P)\SE{Q) = 0. By symmetty, SE{Q)\SE{P) = 0, too. Thus, SE{P) = 

5£;(g). □ 

3.3 Uniform Equivalence and UE-models 

Let D be a set of atoms. By ro we denote a monotone-constraint rule 

rD= {D,{D})^. 

Adding a rule rD to a program forces all atoms in D to be true (independently of the program). 

Monotone-constraint programs P and Q are uniformly equivalent, denoted by P =„ Q, if for 
every set of atoms D, P L) {rn} and Q U {r£,} have the same stable models. 

An SE-model {X, y) of a monotone-constraint program P is a UE-model of P if for every 
SE-model (X', Y) of P with X C X', eitiier X = X' or X' = Y holds. We write UE{P) to 
denote the set of all UE-models of P. Our notion of a UE-model is a generalization of the notion of 
a UE-model due to Liter and Fink (2003) to the setting of monotone-constraint programs. 
Examples. Let us look again at the program we used to illustrate the concept of an SE-model. We 
showed there that (0, {p, q}) and {{p}, {p, q}) are SE-models of P. Directly from the definition of 
UE-models it follows that ({p}, {p, q}) is a UE-model of P. A 

We will now present a characterization of uniform equivalence of monotone-constraint programs 
under the assumption that their sets of atoms are finite. One can prove a characterization of uniform 
equivalence of arbitrary monotone-constraint programs, generahzing one of the results by Eiter and 
Fink (2003). However, both the characterization and its proof are more complex and, for brevity, we 
restrict our attention to the finite case only. 

We start with an auxiliary result, which allows us to focus only on atoms in At{P) when de- 
ciding whether a pair {X, Y) of sets of atoms is an SE-model of a monotone-constraint program 
P. 

Lemma 5. Let P be a monotone-constraint program, X C Y two sets of atoms. Then {X, Y) E 
SE{P) if and only if{X n At{P),Y n At{P)) G SE{P). 

Proof Since X C y is given, and X C y impUes X n At{P) CYD At{P), the first condition of 
the definition of an SE-model holds on both sides of the equivalence. 

Next, we note that for every constraint C, y |= C if and only if y n Dom{C) \= C. Therefore, 
y ^ P if and only \£Yp[At{P) \= P. That is, the second condition of the definition of an SE-model 
holds for (X, Y) if and only if it holds for (X n At{P), Y n At{P)). 

Finally, we observe tiiat P^ = p^nAt(P) p(^x) = P(X n At{P)). Therefore, 

y n hset{p^{X)) = y n hset{p^''^*^^\x n At{P))). 

Since hset{P^^^*^^\X n At{P))) C At{P), it follows tiiat 

Y n hset{P^ {X)) C X 
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if and only if 

Y n At{P) n hset{P^^^^^^\x n At{P))) c x n At{P). 

Thus, X P^ if and only if X n At{P) ^^ynAtiP) P^'^^*(^). That is, the third condition of the 
definition of an SE-model holds for {X, Y) if and only if it holds for (X n At{P), Y n At{P)). □ 

Lemma 6. Let P be a monotone-constraint program such that At{P) is finite. Then for every 
(X, Y) e SE{P) such that X =>^Y, the set 

{X': X CX' C Y, X' ^ Y, {X', Y) e SE{P)} (2) 

has a maximal element. 

Proof. \f At{P)f\X = ^t(P)ny, then for every element y G y \X, is a maximal element 

of the set (2). Indeed, since {X, Y) G SE{P), by Lemma 5, {X n At{P), Y n At{P)) G SE{P). 
Since X n At{P) = F n At{P) and y At{P), X n At{P) = {Y \ [y]) n At{P). Therefore, 
((^\{?/})'~l^t(P),ynAt(P)) G S'E(P). Then from Lemma 5 and the fact C y, we have 

(y \ {y}, Y) G SE{P). Therefore, y \ {y} belongs to the set (2) and so, it is a maximal element 
of this set. 

Thus, let us assume that At{P) n X / At{P) D Y. Let us define X' = X U (y \ At{P)). 
Then X C X' C y and X' 7^ Y. Moreover, no element in X' \ X belongs to At{P). That is, 
X' n At{P) = X n At{P). Thus, by Lemma 5, {X',Y) G SE{P) and so, X' belongs to the set 
(2). Since y \ X' C At{P), by the finiteness of At{P) it follows that the set (2) contains a maximal 
element containing X'. In particular, it contains a maximal element. □ 

Theorem 2. Let P and Q be two monotone-constraint programs such that At{P) U At{Q) is finite. 
Then P =u Q if and only ifUE{P) = UE{Q). 

Proof. (-4=) Let D be an arbitrary set of atoms and y be a stable model of P U {r/j}. Then Y is 
a model of P U {ru}. In particular, y is a model of P and so, {Y, Y) G UE{P). It follows that 
(y, y) G [/ E{Q), too. Thus, y is a model of Q. Since y is a model of r^, P* C y. Consequently, 
y is a model of Q U {ro} and thus, also of {Q U {ro})^ ■ 

Let X = Can{{Q U {ri)})'^, y). Then D Q X CY and, by Proposition 7, X is a y-maximal 
model of {Q U {ro})^- Consequently, X is a y-maximal model of . Since X C y and Y \= Q, 
{X,Y)eSE{Q). 

Let us assume that X ^Y. Then, by Lemma 6, there is a maximal set X' such that X C X' C 

y, X' 7^ y and (X',y) G SE{Q). It follows that (X',y) G [/P(Q). Thus, (X',y) G UE{P) 
and so, X' P^. Since P> C X', X' (P U {rz,})^. We recall that y is a stable model of 
P U {ro}- Thus, y = Can{{P U {r^)})^, y). By Proposition 7, y C X' and so we get X' = Y, 
a contradiction. It follows that X = Y and, consequently, y is a stable model of Q U {vd}- 

By symmetry, every stable model of Q U {rn} is also a stable model of P U {td}- 
(^) First, we note that (Y, Y) G UE(P) if and only if y is a model of P. Next, we note that P and 
Q have the same models. Indeed, the argument used in the proof of Lemma 4 works also under the 
assumption that P Q. Thus, {Y, Y) G UE{P) if and only if {Y, Y) G C/£;(Q). 

Now let us assume that ;7P(P) 7^ UE{Q). Let (X, Y) be an element of {UE{P) \ UE{Q)) U 
(;7P(Q) \ UE{P)). Without loss of generality, we can assume that (X,y) G UE{P) \ UE{Q). 
Since (X, Y) G C/P(P), it follows that 
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1. X C y 

2. Y \= P and, consequently, Y \= Q 

3. X ^Y (otherwise, by our earlier observations, {X, Y) would belong to UE{Q)). 

Let i? = (Q U {rx}V. Clearly, R is a Horn constraint program. Moreover, since Y \= Q and 
X C F, y 1= R. Thus, Can{R, Y) is defined. We have X C Can{R, Y) C Y. We claim that 
Can{R, Y) / Y. Let us assume to the contrary that Can{R, Y) = Y. Then Y G St{Q U {rjc}). 
Hence, Y e St{P U {rx}), that is, y = Can{{P U {rx})'^,!^). By Proposition 7, y is the 
least y-maximal model of {P U {rx})^ and X is a y-maximal model of {P U {^x})^ (since 
{X, Y) G S'^(P), X and so, X (P U {rx}V , too). Consequently, y C X and, as 

X C y, X = y, a contradiction. 

Thus, Can{R, Y) ^ Y. By Proposition 7, Can{R, Y) is a y-maximal model of R. Since 
C R, it follows that Can{R,Y) is a y-maximal model of and so, (Can(i?, y), y) G 
SE{Q). Since Can{R, Y) / y, from Lemma 6 it follows that there is a maximal set X' such that 
CaniR,Y) C X' C Y, X' ^ Y and (X',y) G 5£(Q). By the definition, (X',y) G [/P(g). 
Since (X,y) ^ f/£;(g). X ^ X'. Consequently, since X C X', X' / y and (X,y) G UE{P), 

(x',y) ^ i7£;(P). 

Thus, (X'.y) G UE{Q) \ UE{P). By applying now the same argument as above to {X' ,Y) 
we show the existence of X" such that X' C X" (ZY,X' + X", X" / Y and (X", y) G SE{P). 
Consequently, we have X C X", X ^ X" aadY ^ X", which contradicts the fact that (X, Y) G 
UE{P). It follows then that UE{P) = UE{Q). □ 

Examples. Let P = {l{p,q} not(2{p, g})}, and Q = {p ^ not(g), (7 ^ not(p)}. Then 
P and Q are strongly equivalent. We note that both programs have {p}, {q}, and {p, q} as models. 
Furthermore, {{p}, {p}), {{q}, {q}), {{p}, {p, q}), {{q}, {p, q}), {{p, q}, {p, q}) and (0, {p, q}) are 
"all" SE-models of the two programs ^. 

Thus, by Theorem 1, P and Q are strongly equivalent. 

We also observe that the first five SE-models are precisely UE-models of P and Q. Therefore, 
by Theorem 2, P and Q are also uniformly equivalent. 

It is possible for two monotone-constraint programs to be uniformly but not strongly equivalent. 
If we add rule p ^ to P, and rule p qXo Q, then the two resulting programs, say P' and Q' , are 
uniformly equivalent. However, they are not strongly equivalent. The programs P' U ^ p} and 
Q' U •<— p) have different stable models. Another way to show it is by observing that (0, {p, q}) 
is an SE-model of Q' but not an SE-model of P'. A 

4. Fages Lemma 

In general, supported models and stable models of a logic program (both in the normal case and the 
monotone-constraint case) do not coincide. Fages Lemma (Fages, 1994), later extended by Erdem 
and Lifschitz (2003), establishes a sufficient condition under which a supported model of a nor- 
mal logic program is stable. In this section, we show that Fages Lemma extends to programs with 
monotone constraints. 

3. From Lemma 5 and Theorem 1, it follows that only those SE-models that contain atoms only from At{P) U At{Q) 
are the essential ones. 
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Definition 5. A monotone-constraint program P is called tight on a set M C At{P) of atoms, if 
there exists a mapping Xfrom M to ordinals such that for every rule A ^ y4i, . . . , A^, not(ylfc_|_i), 
. . . ,not(^m) in P(M), if X is the domain of A and Xi the domain of Ai, 1 < i < k, then for 
every x E M f\ X and for every a G M n Uti Xi, \{a) < Xix). 

We will now show that tightness provides a sufficient condition for a supported model to be 
stable. In order to prove a general result, we first establish it in the Horn case. 

Lemma 7. Let P be a Horn monotone-constraint program and let M be a supported model of P. 

If P is tight on M, then M is a stable model of P. 

Proof. Let M be an arbitrary supported model of P such that P is tight on M. Let A be a mapping 
showing the tightness of P on M. We will show that for every ordinal a and for every atom x G M 
such that \{x) < a, x G Can{P, M). We will proceed by induction. 

For the basis of the induction, let us consider an atom x G M such that A(x) = 0. Since M is a 
supported model for P and x G M, there exists a rule r G P{M) such that x G hset{r). Moreover, 
since P is tight on M, for every A G hd{r) and for every y G Dom{A) n M, A(y) < \{x) = 0. 
Thus, for every A G hd{r), Dom{A) n M = 0. Since M \= bd{r) and since P is a Horn monotone- 
constraint program, it follows that |= bd{r). Consequently, hset{r) fl M C Can{P, M) and so, 
X G Can{P,M). 

Let us assume that the assertion holds for every ordinal (5 < a and let us consider x G M such 
that A(x) = a. As before, since M is a supported model of P, there exists a rule r G P{M) such 
that X G hset{r). By the assumption, P is tight on M and, consequently, for every A G bd{r) and 
for every y G Dom{A) n M, X{y) < A(x) = a. By the induction hypothesis, for every A G bd{r), 
Dom{A) n M C Can{P, M). Since P is a Horn monotone-constraint program, Can{P,M) \= 
bd{r). By Proposition 3, hset{r) n M C Can{P, M) and so, x G Can{P, M). 

It follows that M C Can{P, M). By the definition of a canonical computation, we have 
Can{P, M) C M. Thus, M = Can{P, M). By Proposition 6, M is a stable model of P. □ 

Given this lemma, the general result follows easily. 

Theorem 3. Let P be a monotone-constraint program and let M be a supported model of P. If P 
is tight on M, then M is a stable model of P. 

Proof. One can check that if M is a supported model of P, then it is a supported model of the reduct 
P^. Since P is tight on M, the reduct P^ is tight on M, too. Thus, M is a stable model of P^ 
(by Lemma 7) and, consequently, a derivable model of P^ (by Proposition 6). It follows that M is 
a stable model of P. □ 

5. Logic PW^" and the Completion of a Monotone-constraint Program 

The completion of a normal logic program (Clark, 1978) is a propositional theory whose models are 
precisely supported models of the program. Thus, supported models of normal logic programs can 
be computed by means of SAT solvers. Under some conditions, for instance, when the assumptions 
of Pages Lemma hold, supported models are stable. Thus, computing models of the completion 
can yield stable models, an idea implemented in the first version of cmodels software (Babovich & 
Lifschitz, 2002). 
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Our goal is to extend the concept of the completion to programs with monotone constraints. The 
completion, as we define it, retains much of the structure of monotone-constraint rules and allow 
us, in the restricted setting of Iparse programs, to use pseudo-boolean constraint solvers to compute 
supported models of such programs. In this section we define the completion and prove a result 
relating supported models of programs to models of the completion. We discuss extensions of this 
result in the next section and their practical computational applications in Section 8. 

To define the completion, we first introduce an extension of propositional logic with monotone 
constraints, a formalism we denote by PU^^. A formula in the logic PX""^ is an expression built 
from monotone constraints by means of boolean connectives A, V (and their infinitary counterparts), 

and -1. The notion of a model of a constraint, which we discussed earlier, extends in a standard 
way to the class of formulas in the logic PL"^'^. 

For a set L = {^i, . . . , A^, not(Afc+i), . . . , not(^m)} of literals, we define 



Let P be a monotone-constraint program. We form the completion of P, denoted Comp{P), as 
follows: 

1. For every rule r & P v/e include in Comp{P) a PU^'^ formula 



(we note that when the set of rules in P is infinite, the disjunction may be infinitary). 

The following theorem generalizes a fundamental result on the program completion from normal 
logic programming (Clark, 1978) to the case of programs with monotone constraints. 

Theorem 4. Let P be a monotone-constraint program. A set M C Ati^P) is a supported model of 
P if and only if M is a model of Comp{P). 

Proof. (^) Let us suppose that M is a supported model of P. Then M is a model of P, that is, for 
each rule r G P, if M |= bd{r) then M |= hd{r). Since M ^ bd{r) if and only if M |= [bd{r)]^, 
it follows that all formulas in Camp (P) of the first type are satisfied by M. 

Moreover, since M is a supported model of P, M C hset{P{M)). That is, for every atom 
X G M, there exists at least one rule r in P such that x G hset{r) and M \= bd{r). Therefore, all 
formulas in Comp{P) of the second type are satisfied by M, too. 

(<;=) Let us now suppose that M is a model of Comp{P). Since M \= bd{r) if and only if M |= 
[bd{r)]^, and since M satisfies formulas of the first type in Comp{P), M is a model of P. 

Let X G M. Since M satisfies the formula x XJ {[bd{r)\'^ : r E P,x ^ hset{r)}, it follows 
that M satisfies \/ {[bd{r)]^ : r £ P,x G hset{r)}. That is, there is r G P such that M satisfies 
[bd{r)]^ (and so, bd{r), too) and x G hset{r). Thus, x G hset{P{M)). Hence, M is a supported 
model of P. □ 

Theorems 3 and 4 have the following corollary. 



= ^1 A . . . A ^jt A -'Ak+i A ... A -.A 



[M(r)]^ ^ hd{r) 



2. For every atom x G At{P), we include in Comp{P) a PU^'^ formula 
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Corollary 5, Let P be a monotone-constraint program. A set M C. At{P) is a stable model ofP if 
P is tight on M and M is a model of Comp{P). 

We observe that for the material in this section it is not necessary to require that constraints 
appearing in the bodies of program rules be monotone. However, since we are only interested in this 
case, we adopted the monotonicity assumption here, as well. 

6. Loops and Loop Formulas in Monotone-constraint Programs 

The completion alone is not quite satisfactory as it relates supported not stable models of monotone- 
constraint programs with models of PX"^ theories. Loop formulas, proposed by Lin and Zhao 
(2002), provide a way to eliminate those supported models of normal logic programs, which are not 
stable. Thus, they allow us to use SAT solvers to compute stable models of arbitrary normal logic 
programs and not only those, for which supported and stable models coincide. 

We will now extend this idea to monotone-constraint programs. In this section, we will restrict 
our considerations to programs P that aiefinitary, that is, At{P) is finite. This restriction implies 
that monotone constraints that appear in finitary programs have finite domains. 

Let P be a finitary monotone-constraint program. The positive dependency graph of P is the 
directed graph Gp = {V, E), where V = At{P) and (u, v) is an edge in E if there exists a rule 
r £ P such that u G hset{r) and v G Dom(A) for some monotone constraint A G bd{r) (that is, 
A appears non-negated in bd{r)). We note that positive dependency graphs of finitary programs are 
finite. 

Let G = {V, E) be a directed graph. A set L C y is a loop in G if the subgraph of G induced 
by L is strongly connected. A loop is maximal if it is not a proper subset of any other loop in G. 
Thus, maximal loops are vertex sets of strongly connected components of G. A maximal loop is 
terminating if there is no edge in G from L to any other maximal loop. 

These concepts can be extended to the case of programs. By a loop (maximal loop, terminating 
loop) of a monotone-constraint program P, we mean the loop (maximal loop, terminating loop) 
of the positive dependency graph Gp of P. We observe that every finitary monotone-constraint 
program P has a terminating loop, since Gp is finite. 

Let X C At{P). By Gp[X] we denote the subgraph of Gp induced by X. We observe that if 
X 7^ then every loop of Gp[X] is a loop of Gp. 

Let P be a monotone-constraint program P. For every model M of P (in particular, for every 
model M of Comp{P)), we define M~ = M\ Can{P^ , M). Since M is a model of P, M is a 
model of P^ . Thus, Can{P^, M) is well defined and so is M". 

For every loop in the graph Gp we will now define the corresponding loop formula. First, for a 
constraint A = {X, G) and a set L C At, we set A\i = {X, {Y e G : Y n L = 0}) and call A^^ 
the restriction of A to L. Next, let r be a monotone-constraint rule, say 

r= ^ <«- Ai, . . . , Ak,not{Ak+i), not(^^). 

If L C At, then define a PL™'^ formula f^iir) by setting 

/3L(r) = A . . . A Ak\L A -A^+i A ... A ^A^. 

Let L be a loop of a monotone-constraint program P. Then, the loop formula for L, denoted by 
LP{L), is the PX™^ formula 

LP(L) = \/ \/{pLir) -.reP and L n hset{r) / 0} 
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(we recall that we use the convention to write a for the constraint C(a) = ({a}, {{a}}). A loop 
completion of a finitary monotone-constraint program P is the PL'^'^ theory 

LComp{P) = Comp{P) U {LP{L) : L is a loop in Gp}. 

The following theorem exploits the concept of a loop formula to provide a necessary and suffi- 
cient condition for a model being a stable model, transfinite one. 

Theorem 6. Let P be a finitary monotone-constraint program. A set M C At{P) is a stable model 
of P if and only if M is a model of LComp{P). 

Proof. (^) Let M be a stable model of P. Then M is a supported model of P and, by Theorem 4, 

M \= Comp(P). 

Let L be a loop in P. If M n L = then M \= LP{L). Thus, let us assume that M n L / 0. 
Since M is a stable model of P, M is a derivable model of P'^ , that is, M = Can{P'^, M). 
Let {Xn)n=i),i,... be the canonical P^-derivation with respect to M (since we assume that P is 
finite and each constraint in P has a finite domain, P-derivations reach their results in finitely many 
steps). Since Can{P^, M)r\L = MnL^$, there is a smallest index n such that X„ n L / 0. 
In particular, it follows that n > (as Xq = 0) and L D Xn-i = 0. 

Since X„ = hset{P^^ [Xn-i) n M and X„ n L 7^ 0, there is a rule r e P^(X„_i) such that 
hset{r) r\ L ^ that is, such that L fl hset{r)) ^ 0. Let r' be a rule in P, which contributes r to 
P^. Then, for every literal not(A) G 6(f(r'), M ^ not(A). Let .4 G bd{r'). Then ^ € &c?(r) and 
so, Xn-i \= A. Since n L = 0, ^ A\^i^, too. By the monotonicity of A\i^, M \= A^i^. 

Thus, M 1= /3L(r'). Since hset{r') H L ^ $, L n hset{r)) 7^ and so, M ^ LP{L). Thus, 
M ^ LComp{P). 

(<;=) Let us consider a set M C At{P) such that M is not a stable model of P. If M is not a 
supported model of P that M ^ Comp{P) and so M is not a model of LComp{P). Thus, let us 
assume that M is a supported model of P. It follows that M~ ^ 0. Let L C M~ be a terminating 
loopfor Gp[M-]. 

Let r' be an arbitrary rule in P such that L n hset{r')) 7^ 0, and let r be the rule obtained from r' 
by removing negated constraints from its body. Now, let us assume that M \= (5^' (L). It follows that 
for every literal not(^) G bd{r'), M \= not{A). Thus, r G P^. Moreover, since L is a terminating 

loopforGp[M"],forevery constraint A G bd{r'), Dom{A)nM~ C L. Since M |= ^i^, it follows 
that Can{P'^,M) \= A. Consequently, hset{r') n L C hset{r') n M C Can{P'^,M) and so, 
L n Can(P^, M) 7^ 0, a contradiction. Thus, M ^ V{/3r'(-^) : r' G P and L n hset{r')) + 0}. 
Since M ^ V it follows that M ^ LP(L) and so, M ^ LComp{P). □ 

The following result follows directly from the proof of Theorem 6 and provides us with a way 
to filter out specific non-stable supported models from Comp{P). 

Theorem 7. Let P be a finitary monotone-constraint program and M a model of C'omp{P). IfM~ 
is not empty, then M violates the loop formula of every terminating loop ofGp[M~]. 

Finally, we point out that. Theorem 6 does not hold when a program P contains infinitely many 
rules. Here is a counterexample: 
Examples. Let P be the set of following rules: 
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l{ao} <- l{ai} 
l{ai} ^ l{a2} 

Han} ^ l{an+i} 

Let M = {ao, . . . , On, . . .}. Then M is a supported model of P. The only stable model of P is 
0. However, = M \ does not contain any terminating loop. The problem arises because there 
is an infinite simple path in Gp[M~]. Therefore, Gp[M~] does not have a sink, yet it does not have 
a terminating loop either. A 

The results of this section, concerning the program completion and loop formulas — most im- 
portantly, the loop-completion theorem — form the basis of a new software system to compute 
stable models of Iparse programs. We discuss this matter in Section 8. 

7. Programs with Convex Constraints 

We will now discuss programs with convex constraints, which are closely related to programs with 
monotone constraints. Programs with convex constraints are of interest as they do not involve ex- 
plicit occurrences of the default negation operator not, yet are as expressive as programs with 
monotone-constraints. Moreover, they directly subsume an essential fragment of the class of Iparse 
programs (Simons et al., 2002). 

A constraint (X, C) is convex, if for every W,Y,Z <ZX such that C y C Z and M^, Z G C, 
we have y G C. A constraint rule of the form (1) is a convex-constraint rule if ^, ^i, . . . , are 
convex constraints and m = k. Similarly, a constraint program built of convex-constraint rules is a 
convex-constraint program. 

The concept of a model discussed in Section 2 applies to convex-constraint programs. To define 
supported and stable models of convex-constraint programs, we view them as special programs with 
monotone-constraints . 

To this end, we define the upward and downward closures of a constraint A = {X, C) to be 
constraints = (X, C+) and A~ = {X, C~), respectively, where 

C+ = {y C X: for some W eC,W <^Y}, and 
C- = {y C X: for some W ^C,Y <^W]. 

We note that the constraint A'^ is monotone. We call a constraint (X, C) antimonotone if C is closed 
under subset, that is, for every W^, y C X, if y G C and C y then W eC.li is clear that the 
constraint A~ is antimonotone. 

The upward and downward closures allow us to represent any convex constraint as the "con- 
junction" of a monotone constraint and an antimonotone constraint.Namely, we have the following 
property of convex constraints. 

Proposition 8. A constraint {X, C) is convex if and only if C = fl C~. 

Proof (^) Let us assume that C = C+ nC~and let us consider a set M such that M' C M C M" , 
where M',M" G C. it follows that M' G C+ and M" G C". Thus, M G C+ and M G C". 
Consequently, M G C, which implies that {X, C) is convex. 
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(^) The definitions directly imply that C C C+ and C C C". Thus, C C C+ n C". Let us 
consider M G C+ Pi C". Then there are sets M', M" G C such that M' C M and M C M" . Since 
C is convex, M G C. Thus, C+ n C" C C and so, C = C+ n C". □ 

Proposition 8 suggests an encoding of convex-constraint programs as monotone-constraint pro- 
grams. To present it, we need more notation. For a constraint A = {X, C), we call the constraint 
(X, C), where C = V{X) \ C, the dual constraint for A. We denote it by A. It is a direct conse- 
quence of the definitions that a constraint A is monotone if and only if its dual A is antimonotone. 

Let C be a convex constraint. We set mc{C) = {C} if C is monotone. We set ■mc{C) = 
{not(C)}, if C is antimonotone. We define mc{C) = {C"*", not(C~)}, if C is neither monotone 
nor antimonotone. Clearly, C and mc{C) have the same models. 

Let P be a convex-constraint program. By mc{P) we denote the program with monotone con- 
straints obtained by replacing every rule r in P with a rule r' such that 

hd{r') = hd{r)+ and bd{r') = [j{mc{A): A G bd{r)} 

and, if hd{r) is not monotone, also with an additional rule r" such that 

hd{r") = (0,0) and bd{r") = {hd{r)-} U bd{r'). 

By our observation above, all constraints appearing in rules of mc{P) are indeed monotone, that is, 
mc{P) is a program with monotone constraints. 

It follows from Proposition 8 that M is a model of P if and only if M is a model of mc(P). We 
extend this correspondence to supported and stable models of a convex constraint program P and 
the monotone-constraint program mc(P). 

Definition 6. Let P be a convex constraint program. Then a set of atoms M is a supported (or 

stable) model of P if M is a supported (or stable) model of mc(P). 

With these definitions, monotone-constraint programs can be viewed (almost) directly as convex- 
constraint programs. Namely, we note that monotone and antimonotone constraints are convex. 
Next, we observe that if A is a monotone constraint, the expression not(^) has the same meaning 
as the antimonotone constraint A in the sense that for every interpretation M, M \= not (A) if and 
only if M ^ A. 

Let P be a monotone-constraint program. By cc(P) we denote the program obtained from P by 
replacing every rule r of the form (1) in P with r' such that 

hd{r') = hd{r) and bd{r') = [J{Ai: i = 1, . . . ,k} L\[j{A~: j = k + I, . . . ,m} 

One can show that programs P and cc(P) have the same models, supported models and stable 
models. In fact, for every monotone-constraint program P we have P = mc{cc{P)). 
Remark. Another consequence of our discussion is that the default negation operator can be elim- 
inated from the syntax at the price of allowing antimonotone constraints and using antimonotone 
constraints as negated literals. □ 
Due to the correspondences we have established above, one can extend to convex-constraint 
programs all concepts and results we discussed earlier in the context of monotone-constraint pro- 
grams. In many cases, they can also be stated directly in the language of convex-constraints. The 
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most important for us are the notions of the completion and loop formulas, as they lead to new 
algorithms for computing stable models of Iparse programs. Therefore, we will now discuss them 
in some detail. 

As we just mentioned, we could use Comp{mc{P)) as a definition of the completion Comp{P) 
for a convex-constraint logic program P. Under this definition Theorems 9 extends to the case 
of convex-constraint programs. However, Comp{mc{P)) involves monotone constraints and their 
negations and not convex constraints that appear in P. Therefore, we will now propose another 
approach, which preserves convex constraints of P. 

To this end, we first extend the logic PX™'^ with convex constraints. In this extension, which 
we denote by PL'^'^ and refer to as the propositional logic with convex-constraints, formulas are 
boolean combinations of convex constraints. The semantics of such formulas is given by the notion 
of a model obtained by extending over boolean coimectives the concept of a model of a convex 
constraint. 

Thus, the only difference between the logic PL'^^, which we used to define the completion and 
loop completion for monotone-convex programs and the logic PL'^'^ is that the former uses mono- 
tone constraints as building blocks of formulas, whereas the latter is based on convex constraints. In 
fact, since monotone constraints are special convex constraints, the logic PU^*^ is a fragment of the 
logic PL"". 

Let P be a convex-constraint program. The completion of P, denoted by 
Comp{P), is the following set of PL"" formulas: 

1. For every rule r E Pv/e include in Comp{P) a PL"" formula 



(as before, for a set of convex constraints L, denotes the conjunction of the constraints in 



(again, we note that when the set of rules in P is infinite, the disjunction may be infinitary). 

One can now show the following theorem. 

Theorem 8. Let P be a convex-constraint program and let M C At{P). Then M is a supported 

model of P if and only if M is a model of Comp{P). 

Proof. (Sketch) By the definition, M is a supported model of P if and only if M is a supported 
model of mc(P). It is a matter of routine checking that Comp{mc{P)) and Comp{P) have the 
same models. Thus the assertion follows from Theorem 4. □ 

Next, we restrict attention to finitary convex-constraint programs, that is, programs with finite 
set of atoms, and extend to this class of programs the notions of the positive dependency graph and 
loops. To this end, we exploit its representation as a monotone-constraint program mc(P). That is, 
we define the positive dependency graph, loops and loop formulas for P as the positive dependency 
graph, loops and loop formulas of mc(P), respectively. In particular, L is a loop of P if and only if 



[hd{r)f hd{r) 



L) 



2. 



For every atom x G At{P), we include in Comp{P) a PL"" formula 
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L is a loop of mc{P) and the loop formula for L, with respect to a convex-constraint program P, 
is defined as the loop formula LP{L) with respect to the program mc{Py. We note that since loop 
formulas for monotone-constraint programs only modify non-negated Uterals in the bodies of rules 
and leave negated literals intact, there seems to be no simple way to extend the notion of a loop 
formula to the case of a convex-constraint program P without making references to mc{P). 

We now define a loop completion of a finitary convex-constraint program P as the PL"^ theory 

LComp{P) = Comp(P) U {LP{L) : L is a loop of P}. 

We have the following theorem that provides a necessary and sufficient condition for a set of 
atoms to be a stable model of a convex-constraint program. 

Theorem 9. Let P be a finitary convex-constraint program. A set M C. At{P) is a stable model of 
P if and only if M is a model of LComp{P). 

Proof. (Sketch) Since M is a stable model of P if and only of M is a stable model of mc{P), Theo- 
rem 6 impUes that M is a stable model of P if and only if M is a stable model of LComp{mc{P)). 
It is a matter of routine checking that LC omp{mc{P)) and LComp{P) have the same models. 
Thus, the result follows. □ 

In a similar way. Theorem 7 impUes the following result for convex-constraint programs. 

Theorem 10. Let P be a finitary convex-constraint program and M a model of Comp{P). If M" 
is not empty, then M violates the loop formula of every terminating loop of Gp[M~]. 

We emphasize that one could simply use LComp{mc{P)) as a definition of the loop completion 
for a convex-constraint logic program. However, our definition of the completion component of the 
loop completion retains the structure of constraints in a program P, which is important when using 
loop completion for computation of stable models, the topic we address in the next section of the 
paper. 

8. Applications 

In this section, we will use theoretical results on the program completion, loop formulas and loop 
completion of programs with convex constraints to design and implement a new method for com- 
puting stable models of Iparse programs (Simons et al., 2002). 

8.1 Lparse Programs 

Simons et al. (2002) introduced and studied an extension of normal logic programming with weight 
atoms. Formally, a weight atom is an expression 

A = Z[ai = lui, . . . , afe = Wk]u, 

where a^, 1 < z < A; are propositional atoms, and Z, u and Wi, 1 < i < k dXQ non-negative integers. 
If all weights Wi are equal to I, ^ is a cardinality atom, written as l{ai , . . . , a^ju. 

4. There is one minor simplification one might employ. For a monotone constraint A, ^A and A are equivalent and 
A is antimonotone and so, convex. Thus, we can eliminate the operator -i from loop formulas of convex-constraint 
programs by writing A instead of -lA. 
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An Iparse rule is an expression of the form 

A i Ai , . . . , Afi 

where A, Ai, . . . , An are weight atoms. We refer to sets of Iparse rules as Iparse programs. Simons 
et al. (2002) defined for Iparse programs the semantics of stable models. 

A set M of atoms is a model of (or satisfies) a weight atom l[ai = wi, . . . ,ak = Wk]u if 

k 

I < ^{wi : ai e M} < u. 
1=1 

With this semantics a weight atom l[ai = wi, . . . ,ak = Wk]u can be identified with a constraint 
{X, C), where X = {oi, . . . , a^} and 

k 

C = {Y ^X:l< ^{wi : aj e y} < u}. 

i=l 

We notice that all weights in a weight atom W are non-negative. Therefore, if M C M' C M" 
and both M and M" are models of W, then M' is also a model of W. It follows that the constraint 
{X, C) we define above is convex. 

Since (X, C) is convex, weight atoms represent a class of convex constraints and Iparse pro- 
grams syntactically are a class of programs with convex constraints. This relationship extends to 
the stable-model semantics. Namely, Marek and Truszczynski (2004) and Marek et al. (2004, 2006) 
showed that Iparse programs can be encoded as programs with monotone constraints so that the 
concept of a stable model is preserved. The transformation used there coincides with the encoding 
mc described in the previous section, when we restrict the latter to Iparse programs. Thus, we have 
the following theorem. 

Theorem 11. Let P be an Iparse program. A set M Q At is a stable model of P according to the 
definition by Simons et al. (2002) if and only if M is a stable model of P according to the definition 
given in the previous section (when P is viewed as a convex-constraint program). 

It follows that to compute stable models of Iparse programs we can use the results obtained 
earlier in the paper, specifically the results on program completion and loop formulas for convex- 
constraint programs. 

Remark. To be precise, the syntax of Iparse programs is more general. It allows both atoms and 
negated atoms to appear within weight atoms. It also allows weights to be negative. However, 
negative weights in Iparse programs are treated just as a notational convenience. Specifically, an 
expression of the form a = w within a weight atom (where w < G) represents the expression 
not (a) = —w (eliminating negative weights in this way from a weight atom requires modifica- 
tions of the bounds associated with this weight atom). Moreover, by introducing new propositional 
variables one can remove occurrences of negative literals from programs. These transformations pre- 
serve stable models (modulo new atoms). Marek and Truszczynski (2004) and Marek et al. (2004, 
2006) provide a detailed discussion of this transformation. 

In addition to weight atoms, the bodies of Iparse rules may contain propositional literals (atoms 
and negated atoms) as conjuncts. We can replace these propositional literals with weight atoms 
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as follows: an atom a can be replaced with the cardinality atom l{a}, and a literal not(a) — 
with the cardinality atom {a}0. This transformation preserves stable models, too. Moreover, the 
size of the resulting program does not increase more than by a constant factor. Thus, through the 
transformations discussed here, monotone- and convex-constraint programs capture arbitrary Iparse 
programs. □ 

8.2 Computing Stable Models of Lparse Programs 

In this section we present an algorithm for computing stable models of Iparse programs. Our method 
uses the results we obtained in Section 7 to reduce the problem to that of computing models of the 
loop completion of an Iparse program. The loop completion is a formula in the logic PL'^'^, in 
which the class of convex constraints is restricted to weight constraints, as defined in the previous 
subsection. We will denote the fragment of the logic PL'^'^ consisting of such formulas by PL'^"'. 

To make the method practical, we need programs to compute models of theories in the logic 
Piwa ^jjj jjQ^ show a general way to adapt to this task off-the-shelf pseudo-boolean constraint 
solvers (Een & Sorensson, 2003; Aloul et al., 2002; Walser, 1997; Manquinho & Roussel, 2005; Liu 
& Truszczynski, 2003). 

Pseudo-boolean constraints {PB for short) are integer programming constraints in which vari- 
ables have 0-1 domains. We will write them as inequalities 



where comp stands for one of the relations <, >, < and >, Wi's and w are integer coefficients 
(not necessarily non-negative), and XiS are integers taking value or 1. A set of pseudo-boolean 
constraints is a pseudo-boolean theory. 

Pseudo-boolean constraints can be viewed as constraints. The basic idea is to treat each 0-1 vari- 
able X as a prepositional atom (which we will denote by the same letter). Under this correspondence, 
a pseudo-boolean constraint (3) is equivalent to the constraint {X, C), where X = {xi , . . . , x^} and 



in the sense that solutions to (3) correspond to models of (X, C) (xi = 1 in a solution if and only if 
Xi is true in the corresponding model). In particular, if all coefficients wi and the bound w in (3) are 
non-negative, and if comp = '>', then the constraint (3) is equivalent to a monotone lower-bound 
weight atom w[xi = wi, . . . , Xn = Wn]- 

It follows that an arbitrary weight atom can be represented by one or two pseudo-boolean con- 
straints. More generally, an arbitrary PL'^"- formula F can be encoded as a set of PB constraints. 
We will describe the translation as a two-step process. 

The first step consists of converting F to a clausal form Tci{FY. To control the size of the 
translation, we introduce auxiliary propositional atoms. Below, we describe the translation F i-)- 
Tci{F) under the assumption that F is a formula of the loop completion of an Iparse program P. 
Our main motivation is to compute stable models of logic programs and to this end algorithms for 
computing models of loop completions are sufficient. 

5. A PL*"" clause is any formula Bi A . . . A Bm. — >■ -ffi V . . . V Hn, where Bi and Hj are weight atoms. 



wi X xi -\- . . . -\- Wk x Xk comp w, 



(3) 



k 




i=l 
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Let F be a formula in the loop completion of an Iparse-program P. We define Td (F) as follows 
(in the transformation, we use a propositional atom a; as a shorthand for the cardinahty atom C{x) = 
l{x}). 

1. If F is of the form Ai A . . . A ^ A, then TdiF) = F 

2. If F is of the form x {[bd{ri)]^) V ... V {[bd{ri)]^), then we introduce new propositional 
atoms br,i, ... , b^-^i and set Tci{F) to consist of the following PL^"' clauses: 

X br,l V ... V br,l 

[bd{ri)]^ — > hf^i, for every bd{ri) 
br;i — > Aj, for every hd{ri) and Aj € bd{ri) 

3. If F is of the form V L ^ \l r{-I^L{r)}, where L is a set of atoms, and every /^^(r) is a conjunc- 
tion of weight atoms, then we introduce new propositional atoms bdf^^r for every /^^(r) in F and 
represent y L as the weight atom Wl = = 1 : k E L]. We then define Td (F) to consist of the 
following clauses: 

WL^\/bdfL,r 

PUr) bdfL,r, for every ^^^(r) G F 
bdfL,r Aj, for every ^^(r) G F and Aj G /3L{r). 
It is clear that the size Td (F) is linear in the size of F. 

The second step of the translation, converts a PL""* formula (7 in a clausal form into a set of 
PB constraints, Tpii{C). To define the translation C Tpii{C), let us consider a PL^"- clause C of 
the form 

BiA...ABm^Hi\/...\/Hn, (4) 

where i?j's and iJj's are weight atoms. 

We introduce new propositional atoms bi, . . . ,bm and hi, hn to represent each weight atom 
in the clause. As noted earlier in the paper, we simply write x for a weight atoms of the form l[x = 
1]. With the new atoms, the clause (4) becomes a propositional clause b\ A . . . Abm ^ hi\J . . .y hn- 
We represent it by the following PB constraint: 

—bi - . . . — bm + hi + . . . + hn > 1 - m. (5) 

Here and later in the paper, we use the same symbols to denote propositional variables and the cor- 
responding 0-1 integer variables. The context will always imply the correct meaning of the symbols. 
Under this convention, it is easy to see that a propositional clause 6i A . . . A — /ii V . . . V /i„ 
and its PB constraint (5) have the same models. 

We introduce next PB constraints that enforce the equivalence of the newly introduced atoms 
bi (or hi) and the corresponding weight atoms Bi (or Hi). 

Let B = l[ai = wi, . . . ,ak = Wk]u be a weight atom and b a propositional atom. We split B 
to i?+ and B~ and introduce two more atoms b^ and b~. To model B = b,we model with pseudo- 
boolean constraints the following three equivalences: 6 = 6+ A b~, 6+ = B^, and b^ = B^. 
1 . The first equivalence can be captured with three propositional clauses. Hence the following three 
PB constraints model that equivalence: 

-b + b+ >0 (6) 
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-b + b- >0 (7) 
-b+-b~ + b>-l (8) 

2. The second equivalence, 6+ = S+, can be modeled by the following two PB constraints 

k 

(-1) X 6+ + ^(oi xwi)>0 (9) 

i=l 

k k 

-(H Wi-l + l)yib+ + ^{ai xwi)<l-l (10) 

1=1 i=l 

3. Similarly, the third equivalence, b~ = B~, can be modeled by the following two PB constraints 

k k k 

i^Wi - u) X b- + ^{ai X Wi) <^Wi (11) 

1=1 i=l i=l 

k 

{u + l)xb~ + ^(aj xwi)>u + l (12) 
j=i 

We define now Tpb{C), for a PL"'" clause C, as the set of all pseudo-boolean constraints (5) and 
(6), (7), (8), (11), (12), (9), (10) constructed for every weight atom occurring in C. One can verify 
that the size of Tpi,{C) is hnear in the size of C. Therefore, Tpb{Tci{F)) has size hnear in the size of 
F. 

In the special case where all Bi's and Hj's are weight atoms of the form = 1] and l[hj = 1], 
we do not need to introduce any new atoms and PB constraints (6), (7), (8), (11), (12), (9), (10). 
Then Tpb(C) consists of a single PB constraint (5). 

We have the following theorem establishing the correctness of the transformation r. The proof 
of the theorem is straightforward. 

Theorem 12. Let F be a loop completion formula in logic PL^"", and M a set of atoms, M C 
At{F). Then M is a model of F in PU""' logic if and only if M has a unique extension M' by 
some of the new atoms in At(Tpb{Tci{F))) such that M' is a model of the pseudo-boolean theory 

Tpb{Tcl{F)). 

We note that when we use solvers designed for PL™" theories, then translation r^^ is no longer 
needed. The benefit of using such solvers is that we do not need to split weight atoms in the PL'^°- 
theories and do not need the auxiliary atoms introduced in Tpj. 

8.2.1 The Algorithm 

We follow the approach proposed by Lin and Zhao (2002). As in that paper, we first compute the 
completion of a Iparse program. Then, we iteratively compute models of the completion using a 
PB solver. Whenever a model is found, we test it for stability. If the model is not a stable model of 
the program, we extend the completion by loop formulas identified in Theorem 10. Often, adding a 
single loop formula filters out several models of Comp (P) that are not stable models of P. 

The results given in the previous section ensure that our algorithm is correct. We present it in 
Figure 1 . We note that it may happen that in the worst case exponentially many loop formulas are 
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Input: P — a Iparse program; 

A — a pseudo-boolean solver 
BEGIN 

compute the completion Comp{P) of P; 

T:=Tpbirci{Comp{P))y, 

do 

if (solver A finds no models of T) 

output "no stable models found" and terminate; 
M := a model of T found by A; 
if (M is stable) output M and terminate; 
compute the redact P'^ of P with respect to M; 
compute the greatest stable model M', contained in M, of P^; 
M- := M \ M'; 

find all terminating loops in M~; 

compute loop formulas and convert them into PB constraints using 

Tpb and Tci; 

add all PB constraints computed in the previous step to T; 
while (true); 
END 



Figure 1: Algorithm of pbmodels 

needed before the first stable model is found or we determine that no stable models exist (Lin & 
Zhao, 2002). However, that problem arises only rarely in practical situations^. 

The implementation of pbmodels supports several PB solvers such as satzoo (Een & Sorensson, 
2003), pbs (Aloul et al, 2002), wsatoip (Walser, 1997). It also supports a program wsatcc (Liu & 
Truszczynski, 2003) for computing models of PL'""^ theories. When this last program is used, the 
transformation, from "clausal" PL""* theories to pseudo-boolean theories is not needed. The first 
two of these four programs are complete PB solvers. The latter two are local-search solvers based 
on wsat (Selman, Kautz, & Cohen, 1994). 

We output the message "no stable model found" in the first line of the loop and not simply "no 
stable models exist" since in the case when A is a local-search algorithm, failure to find a model 
of the completion (extended with loop formulas in iteration two and the subsequent ones) does not 
imply that no models exist. 

8.3 Performance 

In this section, we present experimental results concerning the performance of pbmodels. The ex- 
periments compared pbmodels, combined with several PB solvers, to smodels (Simons et al, 2002) 
and cmodels (Babovich & Lifschitz, 2002). We focused our experiments on problems whose state- 



6. In fact, in many cases programs turn out to be tight with respect to their supported models. Therefore, supported 
models are stable and no loop formulas are necessary at all. 
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ments explicitly involve pseudo-boolean constraints, as we designed pbmodels with such problems 
in mind. 

For most benchmark problems we tried cmodels did not perform well. Only in one case (vertex- 
cover benchmark) the performance of cmodels was competitive, although even in this case it was 
not the best performer. Therefore, we do not report here results we compiled for cmodels. For a 
complete set of results we obtained in the experiments we refer to http : / / www .cs.uky.edu/ 

ai/pbmodels. 

In the experiments we used instances of the following problems: traveling salesperson, weighted 
n-queens, weighted Latin square, magic square, vertex cover, and Towers of Hanoi. The Iparse 
programs we used for the first four problems involve general pseudo-boolean constraints. Programs 
modeling the last two problems contain cardinality constraints only. 

Traveling salesperson problem (TSP). An instance consists of a weighted complete graph with n 
vertices, and a bound w. All edge weights and w are non-negative integers. A solution to an instance 
is a Hamiltonian cycle whose total weight (the sum of the weights of all its edges) is less than or 
equal to w. 

We randomly generated 50 weighted complete graphs with 20 vertices, To this end, in each case 
we assign to every edge of a complete undirected graph an integer weight selected uniformly at 
random from the range [1..19]. By setting w to 100 we obtained a set of "easy" instances, denoted 
by TSP-e (the bound is high enough for every instance in the set to have a solution). From the same 
collection of graphs, we also created a set of "hard" instances, denoted by TSP-h, by setting w to 62. 
Since the requirement on the total weight is stronger, the instances in this set in general take more 
time. 

Weighted n-queens problem (WNQ). An instance to the problem consists of a weighted n x n 
chess board and a bound w. All weights and the bound are non-negative integers. A solution to an 
instance is a placement of n queens on the chess board so that no two queens attack each other and 
the weight of the placement (the sum of the weights of the squares with queens) is not greater than 

w. 

We randomly generated 50 weighted chess boards of the size 20 x 20, where each chess board 
is represented by a set of n x n integer weights Wij, 1 < i,j < n, all selected uniformly at random 
from the range [1..19]. We then created two sets of instances, easy (denoted by wnq-e) and hard 
(denoted by wnq-h), by setting the bound w to 70 and 50, respectively. 

Weighted Latin square problem (WLSQ). An instance consists of an n x n array of weights Wi^j, 
and a bound w. All weights Wi^j and w are non-negative integers. A solution to an instance is an 
n X n array L with all entries from {1, . . . , n} and such that each element in {1, . . . , n} occurs 
exactly once in each row and in each column of L, and Y17=i '^"i=i ^I^'i] ^ "^ij — 

We set ri = 10 and we randomly generated 50 sets of integer weights, selecting them uniformly 
at random from the range [1..9]. Again we created two famihes of instances, easy {wlsq-e) and hard 
(wlsq-h), by setting w to 280 and 225, respectively. 

Magic square problem. An instance consists of a positive integer n. The goal is to construct an 
n X n array using each integer 1, ... as an entry in the array exactly once in such a way that 
entries in each row, each column and in both main diagonals sum up to n{n? + l)/2. For the 
experiments we used the magic square problem for n = 4,5 and 6. 

Vertex cover problem. An instance consists of graph with n vertices and m edges, and a non- 
negative integer k — a bound. A solution to the instance is a subset of vertices of the graph with no 
more than k vertices and such that at least one end vertex of every edge in the graph is in the subset. 
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We randomly generated 50 graphs, each with 80 vertices and 400 edges. For each graph, we set 
A; to be a smallest integer such that a vertex cover with that many elements still exists. 

Towers of Hanoi problem. This is a slight generalization of the original problem. We considered 
the case with six disks and three pegs. An instance consists of an initial configuration of disks that 
satisfies the constraint of the problem (larger disk must not be on top of a smaller one) but does not 
necessarily requires that all disks are on one peg. These initial configurations were selected so that 
they were 31, 36, 41 and 63 steps away from the goal configuration (all disks from the largest to the 
smallest on the third peg), respectively. We also considered a standard version of the problem with 
seven disks, in which the initial configuration is 127 steps away from the goal. 

We encoded each of these problems as a program in the general syntax of Iparse, which allows 
the use of relation symbols and variables (Syrjanen, 1999). The programs are available at http : 
//www. cs . uky . edu/ai/pbmodels. We then used these programs in combination with ap- 
propriate instances as inputs to Iparse (Syrjanen, 1999). In this way, for each problem and each 
set of instances we generated a family of ground (prepositional) Iparse programs so that stable 
models of each of these programs represent solutions to the corresponding instances of the prob- 
lem (if there are no stable models, there are no solutions). We used these families of Iparse pro- 
grams as inputs to solvers we were testing. All these ground programs are also available at http : 
//www . cs . uky . edu/ai/pbmodels. 

In the tests, we used phmodels with the following four PB solvers: satzoo (Een & Sorensson, 
2003), pbs (Aloul et al., 2002), wsatcc (Liu & Truszczynski, 2003), and wsatoip (Walser, 1997). In 
particular, wsatcc deals with PL'"'"' theories directly. 

All experiments were run on machines with 3.2GHz Pentium 4 CPU, 1GB memory, running 
Linux with kernel version 2.6.11, gcc version 3.3.4. In all cases, we used 1000 seconds as the 
timeout limit. 

We first show the results for the magic square and towers of Hanoi problems. In Table 1, for 
each solver and each instance, we report the corresponding running time in seconds. Local-search 
solvers were unable to solve any of the instances in the two problems and so are not included in the 
table. 



Benchmark 


smodels 


pbmodels-satzoo 


pbmodels-pbs 


magic square (4 x 4) 


1.36 


1.70 


2.41 


magic square (5 x 5) 


> 1000 


28.13 


0.31 


magic square (6 x 6) 


> 1000 


75.58 


> 1000 


towers of Hanoi {d — 6,t = 31) 


16.19 


18.47 


1.44 


towers of Hanoi {d = 6,t = 36) 


32.21 


31.72 


1.54 


towers of Hanoi {d = 6,t = 41) 


296.32 


49.90 


3.12 


towers of Hanoi (d — 6,t = 63) 


> 1000 


> 1000 


3.67 


towers of Hanoi {d = 7,t = 127) 


> 1000 


> 1000 


22.83 



Table 1 : Magic square and towers of Hanoi problems 

Both pbmodels-satzoo and pbmodels-pbs perform better than smodels on programs obtained 
from the instances of both problems. We observe that pbmodels-pbs performs exceptionally well 
in the tower of Hanoi problem. It is the only solver that can compute a plan for 7 disks, which 
requires 127 steps. Magic square and Towers of Hanoi problems are highly regular. Such problems 
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# of SAT instances 


# of UNSAT instances 


# of UNKNOWN instances 


TSP-e 


50 








TSP-h 


31 


1 


18 


wnq-e 


49 





1 


wnq-h 


29 





21 


wlsq-e 


45 


4 


1 


wlsq-h 


8 


41 


1 


vtxcov 


50 









Table 2: Summary of Instances 





smodels 


pbmodels-satzoo 


pbmodels-pbs 


TSP-e 


45/17 


50/30 


18/3 


TSP-h 


7/3 


16/14 


0/0 


wnq-e 


11/5 


26/23 


0/0 


wnq-h 


2/2 


0/0 


0/0 


wlsq-e 


21/1 


49/29 


46/19 


wlsq-h 


0/0 


47/26 


47/23 


Ytxcoy 


50/ 10 


50/1 


17/:5 


sum over all 


136/68 


238/123 


158/48 



Table 3: Summary on all instances 



are often a challenge for local-search problems, which may explain a poor performance we observed 
for pbmodels-wsatcc and pbmodels-wsatoip on these two benchmarks. 

For the remaining four problems, we used 50-element families of instances, which we gen- 
erated randomly in the way discussed above. We studied the performance of complete solvers 
(smodels, pbmodels-satzoo and pbmodels-pbs) on all instances. We then included local-search 
solvers (pbmodels-wsatcc and pbmodelswsatoip) in the comparisons but restricted attention only 
to instances that were determined to be satisfiable (as local-search solvers are, by their design, un- 
able to decide unsatisfiability). In Table 2, for each family we list how many of its instances are 
satisfiable, unsatisfiable, and for how many of the instances none of the solvers we tried was able to 
decide satisfiability. 

In Table 3, for each of the seven families of instances and for each complete solver, we report 
two values s/w, where s is the number of instances solved by the solver and w is the number of 
times it was the fastest among the three. 

The results in Table 3 show that overall pbmodels-satzoo solved more instances than pbmodels- 
pbs, followed by smodels. When we look at the number of times a solver was the fastest one, 
pbmodels-satzoo was a clear winner overall, followed by smodels and then by pbmodels-pbs. 
Looking at the seven families of tests individually, we see that pbmodels-satzoo performed better 
than the other two solvers on five of the famiUes. On the other two smodels was the best performer 
(although, it is a clear winner only on the vertex-cover benchmark; all solvers were essentially 
ineffective on the wnq-h). 

We also studied the performance of pbmodels combined with local-search solvers wsatcc (Liu 
& Truszczyriski, 2003) and wsatoip (Walser, 1997). For this study, we considered only those in- 
stances in the seven families that we knew were satisfiable. Table 4 presents results for all solvers 
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smodels 


pbmd-satzoo 


pbmd-pbs 


pbmd-wsatcc 


pbmd-wsatoip 


TSP-e 


45/3 


50/5 


18/2 


32/7 


47/34 


TSP-h 


7/0 


16/2 


0/0 


19/6 


28/22 


wnq-e 


11/0 


26/0 


0/0 


49/45 


49/4 


wnq-h 


2/0 


0/0 


0/0 


29/15 


29/14 


wlsq-e 


21/0 


45/0 


44/0 


45/33 


45/14 


wlsq-h 


0/0 


7/0 


8/0 


7/1 


8/7 


vtxcov 


50/0 


50/0 


47/0 


50/36 


50/15 


sum aver all 


136/3 


194/7 


117/2 


231/143 


256/110 



Table 4: Summary on SAT instances 



we studied (including the complete ones). As before, each entry provides a pair of numbers s/w, 
where s is the number of solved instances and w is the number of times the solver performed better 
than its competitors. 

The results show superior performance of pbmodels combined with local-search solvers. They 
solve more instances than complete solvers (including smodels). In addition, they are significantly 
faster, winning much more frequently than complete solvers do (complete solvers were faster only 
on 12 instances, while local-search solvers were faster on 253 instances). 

Our results demonstrate that pbmodels with solvers of pseudo-boolean constraints outperforms 
smodels on several types of search problems involving pseudo-boolean (weight) constraints). 

We note that we also analyzed the run-time distributions for each of these families of instances. 
A run-time distribution is regarded as a more accurate and detailed measure of the performance of 
algorithms on randomly generated instances^. The results are consistent with the summary results 
presented above and confirm our conclusions. As the discussion of run-time distributions requires 
much space, we do not include this analysis here. They are available at the website http : / / www . 
cs . uky . edu/ai/pbmodels. 

9. Related work 

Extensions of logic programming with means to model properties of sets (typically consisting of 
ground terms) have been extensively studied. Usually, these extensions are referred to by the com- 
mon term of logic programming with aggregates. The term comes from the fact that most properties 
of sets of practical interest are defined through "aggregate" operations such as sum, count, maxi- 
mum, minimum and average. We chose the term constraint to stress that we speak about abstract 
properties that define constraints on truth assignments (which we view as sets of atoms). 

Mumick, Pirahesh, and Ramakrishnan (1990), and Kemp and Stuckey (1991) were among the 
first to study logic programs with aggregates. Recently, Niemela et al. (1999) and Simons et al. 
(2002) introduced the class of Iparse programs. We discussed this formalism in detail earlier in this 
paper. 

Pelov (2004) and Pelov et al. (2006) studied a more general class of aggregates and devel- 
oped a systematic theory of aggregates in logic programming based on the approximation theory 
(Denecker, Marek, & Truszczynski, 2000). The resulting theory covers not only the stable models 
semantics but also the supported-model semantics and extensions of 3-valued Kripke-Kleene and 

7. Hoos and Stiitzle (2005) provide a detailed discussion of this matter in the context of local-search methods. 
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well-founded semantics. The formalism introduced and studied by Pelov (2004) and Pelov et al. 
(2006) allows for arbitrary aggregates (not only monotone ones) to appear in the bodies of rules. 
However, it does not allow for aggregates to appear in the heads of program clauses. Due to dif- 
ferences in the syntax and the scope of semantics studied there is no simple way to relate Pelov' s 
(2004) and Pelov et al.'s (2006) formalism to programs with monotone (convex) constraints. We 
note though that programs with abstract monotone constraints with the heads of rules of the form 
C (a) can be viewed almost literally as programs in the formalism by Pelov (2004) and Pelov et al. 
(2006) and that they have the same stable models according to the definitions we used in this paper 
and those by Pelov (2004) and Pelov et al. (2006). 

Faber et al. (2004) developed the theory of disjunctive logic programs with aggregates. Similarly 
as Pelov (2004) and Pelov et al. (2006), Faber et al. (2004) do not allow for aggregates to appear 
in the heads of program clauses. This is one of the differences between that approach and programs 
with monotone (convex) constraints we studied here. The other major difference is related to the 
postulate of the minimality of stable models (called answer sets in the context of the formalism 
considered by Faber et al., 2004). In keeping with the spirit of the original answer-set semantics 
(Gelfond & Lifschitz, 1991), answer sets of disjunctive programs with aggregates, as defined by 
Faber et al. (2004), are minimal models. Stable models of programs with abstract constraints do 
not have this property. However, for the class of programs with abstract monotone constraints with 
the heads of rules of the form C{a) the semantics of answer sets defined by Faber et al. (2004) 
coincides with the semantics of stable models by Marek and Truszczyhski (2004) and Marek et al. 
(2004, 2006). 

Yet another approach to aggregates in logic programming was presented by Son and PontelU 
(2006). That approach considered programs of the syntax similar to programs with monotone ab- 
stract constraints. It allowed arbitrary constraints (not only monotone ones) but not under the scope 
of not operator. A general principle behind the definition of the stable-model semantics by Son 
and Pontelli (2006) is to view a program with constraints as a concise representation of a set of its 
"instances", each being a normal logic program. Stable models of the program with constraints are 
defined as stable models of its instances and is quite different from the operator-based definition 
by Marek and Truszczyhski (2004) and Marek et al. (2004, 2006). However, for programs with 
monotone constraint atoms which fall in the scope of the formaUsm of Son and PontelU (2006) both 
approaches coincide. 

We also note that recently Son et al. (2006) presented a conservative extension of the syntax 
proposed by Marek and Truszczyhski (2004) Marek et al. (2006), in which clauses are built of 
arbitrary constraint atoms. 

Finally, we point out the work by Ferraris and Lifschitz (2004) and Ferraris (2005) which treats 
aggregates as nested expressions. In particular, Ferraris (2005) introduces a propositional logic with 
a certain nonclassical semantics, and shows that it extends several approaches to programs with 
aggregates, including those by Simons et al. (2002) (restricted to core Iparse programs) and Faber 
et al. (2004). The nature of the relationship of the formaUsm by Ferraris (2005) and programs with 
abstract constraints remains an open problem. 

10. Conclusions 

Our work shows that concepts, techniques and results from normal logic programming, concerning 
strong and uniform equivalence, tightness and Fages lemma, program completion and loop formu- 
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las, generalize to the abstract setting of programs with monotone and convex constraints. These 
general properties specialize to new results about Iparse programs (with the exception of the char- 
acterization strong equivalence of Iparse programs, which was first obtained by Turner, 2003). 

Given these results we implemented a new software pbmodels for computing stable models of 
Iparse programs. The approach reduces the problem to that of computing models of theories con- 
sisting of pseudo-boolean constraints, for which several fast solvers exist (Manquinho & Roussel, 
2005). Our experimental results show that pbmodels with PB solvers, especially local search PB 
solvers, performs better than smodels on several types of search problems we tested. Moreover, as 
new and more efficient solvers of pseudo-boolean constraints become available (the problem is re- 
ceiving much attention in the satisfiability and integer programming communities), the performance 
of pbmodels will improve accordingly. 
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